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Abstract.  By  a  method  which  modifies  Dingle's  original  one, 
broadening  effects  on  the  magnvLuthermal  oscillations  in  two- 
dimensional  electron  systems  are  created  analytically  and 
relevant  thermodynamic  quantities  are  evaluated  explicitly. 

In  particular,  the  behaviors  of  the  cJjcmlcal  potential  and  specific 
heat  in  a  strong  magnetic  field  are  investigated.  Landau 
level  broadening  does  not  affect  the  period  ■  but  causes 
significant  reduction  of  the  amplitude  of  the  oscillations. 

For  its  sensitivity,  the  oscillating  pattern  is  a  good  indicator 
of  level  broadening,  while  the  period  of  oscillation  can  be 
used  for  the  determination  of  the  effective  Bohr  magneton. 
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Introduction 


Two-dimensional  electron  systems  which  are  formed  in  Si  inversion  layers  or 
GaAs/GaAlAs  superlattices  show  very  unusual  properties  especially  In  strong 
magnetic  field.  While  these  properties  have  mostly  been  studied  under  the 
Isothermal  conditions.  It  has  been  pointed  out  (Islhara  and  Shlwa  1983 
Zawadskl  and  Lassni/^g  1984)  that  their  temperature  oscillates  very  strongly 
when  the  magnetic  field  is  varied  adiabatlcally.  The  significance  of 
these  oscillations  Is  due  to  the  lack  of  electron  notion  In  the  direction  of 
the  iMgnetlc  field.  This  makes  the  conversion  of  the  field  energy  Into 
the  kinetic  energy  of  the  electrons  effective. 

Zn  our  previous  work,  hereafter  to  be  called  1  (Islhara  and  Shlwa  1984) 
we  Investigated  the  case  without  level  broadening.  Since  In  actual  systems. 
Impurity  scatterings  cause  level  broadening,  we  give  In  the  present  article 
a  somewhat  comprehensive  treatment  of  the  case  with  level  broadening. 
Different  from  Zawadzkl  and  Lassnig  who  employed  a  numerical  approach,  we 
shall  make  an  analytical  approach  to  the  magnetothermal  ef'  ':t  and  derive 
explicit  formulae  for  relevant  physical  quant  it ites.  ■  For  this  purpose, 
we  shall  employ  a  method  which  modifies  Dingle's  (1952)  for  the  de  Haas- 
van  Alphen  effect.  'This  is  a  phenomenological  approach,  but  gives  the 
advantage  of  being  applicable  to  several  different  cases.  Therefore,  we 
shall  treat  the  cases  of  Lorentzlan  as  well  as  elliptic  broadenlngs.  For 
2D  electron  systems,  Che  latter  type  of  broadening  has  been  Introduced 
effectively  by  Ando  and  Uemura  (1984).  Although  the  final  broadening 
effects  depend  on  the  iMgnlCude  of  the  respective  broadening  parameter. 

It  Is  generally  considered  that  elliptic  broadening  describes  low  temperature 
phenomena  well. 

'  In  actual  2D  ayatems,  there  are  Coulomb  Interactions.  We  remark  chat 
Islhara,  Tsai  and  Wadatl  (1971)  showed  that  these  Interactions  lead  to 
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«  reduction  of  the  amplitude  of  3D  dHvA  oscillations  similarly  to  the  case  of 
Dingle.  More  recenlcy,  Shi«fa  and  Islhara  (1983)  treated  dHvA  oscillations 
in  2D  electron  systems  with  Coulomb  Interaction. 

^  The  magnetothermal  effect  requires  the  evaluation  of  the  entropy  at 
finite  temperature.  This  is  a  difficult  task  In  the  presence  of  Impurity 
scattering  and  a  strong  magnetic  field.  However,  the  Dingle  method  can  be 
used  effectively. 

The  energy  levels  of  2D  electrons  can  be  written  as 

c^-  (2n+l)UgH  f  •fi'/H,  n  -  0,  1,  2,  ...  (1.1) 

where  ii_  >  e)l/2mc  Is  the  effective  Bohr  magneton  with  the  effective  mass  m, 

O 

■  eK/2n  c  is  the  real  Bohr  magneton  with  the  bare  electron  mass  m  ,  and 
B  .  o  o 

g  is  the  effective  Lande's  g  factor,  H  being  the  magnetic  field. 

We  note  that  the  magnetothermal  oscillations  are  primarily  due  to  the 
orbital  motion  of  the  electrons  as  in  the  case  of  the  dHvA  effect.  Therefore, 

we  shall  start  with  the  case  without  Che  spin-magnetic  field  coupling  energy. 

P 

As  can  be  guessed,  this  coullng  causes  a  shift  in  the  phase  of  the  oscillations. 

A 

In  Che  next  section,  we  shall  derive  a  new  formula  for  the  grand 

potential  of  a  20  electron  system  with  level  broadening  by  a  mehtod  which 

modifies  slightly  the  original  Dingle's.  Section  3  gives  a  basic  formula  for 

the  SMgnetoCherawl  effect  for  arbitrary  broadening  and  low  but  finite  temperatuers. 

« 

In  Section  4,  we  present  an  explicit  limiting  expression  for  the  'sMigneCoChermal 

* 

effect  at  absolute  zero.  At  the  saise  time,  we  proscMit  a  new  specific  heat 
formula  for  Che  case  with  level  broadening.  As  we  shall  see,  Che  behavior 
of  the  specific  heat  Is  very  crucial  to  the  magnetothermal  effect.  Finally, 
in  Section  5  we  present  explicit  numerical  results  and  discussions. 

V 
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2.  Effect*  of  Level  Broadening  on  the  Grand  Potential 

In  this  section*  we  treat  level  broadening  by  a  method  which  laodlfies 
slightly  Dingle's  original  theory.  For  analytical  simplicity,  let  us 
consider  an  idealized  2D  case  in  which  Is  set  equal  to  so  that  the 


grand  pot^^ial  becomes 

ft(li)  -  -kT  X«n(l+ 


(2.1) 


where  B  ■  1/kT  is  the  reciprocal  thermal  energy  and  g  Is  the  chemical 


potential.  We  write 


where 


fi(g)  -  I  n^(g+ggH)  +  ifl^Cg-ggH) 


ft  (g)  -  -kT  I  ln[l  + 
n 


(2.2) 


(2.3) 


c”  -  (2ntl)g_H 
n  B 


(2.4) 


Since  the  two  terms  In  equation  (2.2)  differ  from  each  other  only  in  the 


sign  of  y.H,  we  derive  a  formula  for  ft  (p).  This  is  then  equivalent  to 

0  o 

treating  the  case  without  the  spin  contribution. 

We  Introduce  the  density  of  states  #g(<^)  Par  unit  energy  Interval 


given  by 


4  (c)  -  2g  c®) 


where  A 'is  the  surface  area,  re  being  the  electron  mass.  We  express 


equation  (2.3)  by  an  Integral 


ft^(g)  -  -kT  I  de  ♦^(e)  tn(l+e 

Integrating  twice  by  parts,  we  find 

3f  (r) 

fto(w)  -  I  dt  ♦,(01  -jf—  J 


(<(ii-i;) , 


(2.5) 


(2.6) 


(2.7) 
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t^erc 


♦^(e)  -  J  de'l  de" 


o  o 


‘o«>  ■  e(c'-i!)— 

e  +1 


(2.9) 


So  far  the  calculation  Is  formally  exact.  We  now  replace  Dirac's 
delta  functions  by  Lorentzian  functions  of  a  half  width  •ti/T: 


...  AO  »  1  _ 1/T _ 

°  ^  ^  n-0  ^  (c-c  °)^+(l/T)^ 


(2.10) 


where,  and  what  follows  in  this  article,  we  have  setdT-  1  and  2m  -1, 
and  a  -  ^3^* 

At  low  temperatures,  falls  off  very  rapidly  on  either  sides 

of  y.  For  6y»l  and  WT»1,  the  n  sum  may  Include  all  negative  integers, 
enabling  the  use  of  the  Poisson  summation  formula# 

The  details  of  the  calculation  and  a  comparison  with  the  original  Dingle's 

ccJ 

method  are  given  respectively  in  Appendix  A  and  B. 

The  general  case  in  which  the  electron  spins  are  Included  can  be 
treated  in  a  similar  way.  '  The  corresponding  grand  potential  is  given 
for  arbitrary  g  by 


jo'  ^  ^ 


(2.11) 


m  fxi  t  co8(Tri/Y)  C08(3irt) 

+  ^  t  (-)‘^^wSr)  - = - ?— } 

,  n  1-1  I  8inh(iT^£/a) 

u) J  =  P/'  .  _  __  .  _ . 

As  we  see,  the  broadening  effects  appear  only  through  the  functions 
D(r)  and  W(r)  which  are  defined  by 


iHij 


.'»v  -.i'.  v 


V 


!7T!r?T5TrrtT77'="7Tr^TJT^ 


W(r)  -  exp[-iT/(Ta  )  ) 


o,r,  -  if  " 

TI  k=l  k 


2Ta 


(2.12) 


H  (2.13) 


Note  that  W(r)  Is  simply  the  Fourier  transform  of  the  characteristic 
Lorentzlan  form  In  the  sum  of  equation  (2.10).  As  each  energy  level  Is 
broadened  In  accordance  with  this  form  In  energy  space*  the  contribution 
from  the  level  to  the  grand  potential  is  reduced  by  the  amplitude,  reduction 
factor  W(r)  such  that  the  higher  harmonic  contributions  are  progressively 
reduced.  Roreover,  level  broadening  affects  the  non-oscillating  part 
of  the  grand  potential  through  the  function  D(r).  Note  that 


n(0) 


(2.14) 


so  that  the  non-oscillating  part  is  reduced  to  that  given  by  Islhara  and 

f 

Kojima  (1979).  As  broadening  increases*  that  is  when  F  Incrases*  D(r) 

A 

decreases.  This  means  that  the  spin  paramagnetic  contribution  becomes 
relatively  more  effective,  while  the  oscillating  contribution  due  to  the 
orbital  motion  of  the  electrons  is  weakened. 
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3.  Impurity  Effect*  on  Magnetothermal  Effect 

The  temperature  change*  under  adiabatic  variations  of  magnetic 
field  are  given  by 

,3. 

(^) 

'3t  'H.n 

where  n  is  the  electron  density, and  the  entropy  S  is  related  to  the  grand 
partition  function  by 


S .  |j 


(3.2) 


•  ■  ■  (^  fi) 


The  grand  potential  12  ha*  been  evaluated  as  in  equation  (2.29).  In  ^rder 
to  perform  the  differentiation*  of  the  entropy  as  in  equation  (3.1)»  it  1* 
necessary  to  know  the'  behavior  of  the  relaxation  time  t  which  ha*  been 
introduced  phenomenologically.  Let  us  assume  that  t  is  independent  of  field 
and  temperature,  although  it  may  depend  on  the  electron  and  impurity 
concentre  t ion* . 

Since  the  electron  density  is  kept  constant  in  the  derivatives  in  equai- 
tion  (3.1) (While  the  grand  potential  has  been  obtained  for  a' given  chemical 
potential  U,  we  need  the  relation  for  the  number  density  n: 


(3.3) 


T^ils  yields  the  inqxjrtant  relation  between  Y  and 


1.  .  i  .  211  y 
Yo  y  “  t-i 


-  8in{^)cos(|irJl) 
_ —  T  * 


M  (D 


sinh  ( - ) 

a 


(3.4) 


'  ^  < 


where  »  1,  2m  ■  1) 


» 


1_ 

Y, 


2Vn 

2 


1 

Y 


are  dimensionless  parameters.  Note 

variable  such  as  n  or  H,  while  1/y 
2 

the  field  energy  a  . 

For  g  -  2  and  T  •=  0,  equation 

1  1  2  _  w*'  , 

-  .  ~  -  I  ^  sin(-) 


(3.5) 

that  1/Yq  represents  the  experimental 
measures  the  chemical  potential  by 

(3. A)  becomes 


(3.6) 


1 

Y 


2  -Ir 
+  -  tan  [ 

tT 


w  sin(^) 
1-W  cos(^) 


■] 


The  latter  formula  is  valid  for  W  «  exp(-2TTr)<l. 

For  g  «  2,  the  entropy  is  obtained  in  the  following  form 


(3.7) 


(3.8) 


4TtS 

Ak 


2TI 

38 


where 


(3.9) 


L  (x)  «  cothx  -  i  (3.10) 

is  the  Langevin  function.  It  is  worth  noting  that  both  l/Y^  of  equation  (3. A) 
and  S  of  equation  (3.9)  do  not  have  the  term  D(r )  which  appears  In  the  grand 
potential  and  that  the  effects  of  level  broadening  appear  only  through 
the  function  W(r). 

We  conpute  now  the  two  relevant  derivatives  of  the  entropy. 


iMt  UB  introduce 


_ _ !— 

H  <3H  ’T.n 


T  2  ^'3T  'H,n 

Aka 


After  straightforward  but  somewhat  lengthy  calculations,  which  require  the 

!  • 

derivatives  of  the  chemical  potential,  we  arrive  at 


S„  -  r  ■Al'^'(l')-(l'cothl'+int'L(i')] 
H  ^  I 


-C2trs,(l'cothi'+£f)+C]Cl+2lA'^,3*^‘£s/'L(i') 
1  *■  A  *■  A 


T  6  J  A 


-25CE8-  A'  L(  A')  ]^*  [1+2EA"cb 
A  *■  A 

where  for  simplicity  we  heve  used 


f  r  I  t'  9F 


-  A  cos  (AC)  "I  _  sin  (Ac)  f  ,  ^_P  _  ^ 

®A  -  s’inlTI^  *  'a  "  siniTT-  ’  ^  -  27rr  - 

The  Buignctotheraul  effect  is  expressed  in  a  dimensionless  way  by 


»  fli  M 

T  "  S_  H 
T 

FoziBulae  (3.13)  and  (3.14)  are  derived  for  low  Init  finite  temperatures 


The  formula  in  equation  (3.14)  can  be  used  to  derive  the  electronic 


specific  heat  at  constant  area  and  magnetic  field  through 

2 


ka 


'A,H 


(3.17) 


The  first  term  in  equation  (3.14)  yields  the  linear  specific  heat 

2 


c°  .^T 
A  6 


(3.18) 


which  is  the  correct  limiting  expression  in  the  absence  of  a  magnetic 
field. 


In  an  arbitrary  magnetic  field,  the  specific  heat  decreases 
exponentially  with  temperature  (Shlwa  a.id  Islhara  1984).  However, 
so  long  as  r ,  the  broadening  parameter  defined  by  equation  (2.13),  does 


not  vanish,  the  linear  specific  heat  is  retained.  Indeed,  we  find 

(  <  C  ) 

2 


near  absolute  zero 


'A,H 


— T[  1  +  2W  ~ 


1-2Wco8(it/y)h 


(3.19) 
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Magnetothermal  Effect  at  Absolute  Zero 

The  magnetothermal  effect  as  expressed  by  the  dimensionless 


quantity  HdT/ (TdH)  depends  on  the  two  entropy  derivatives  S„  and  S 

H  i 

given  by  equations  (3.11)  and  (3.12)  respectively.  Note  that  these 

two  quantities  are  also  dimensionless.  In  the  natural  units  which  have 

2 

been  adopted  In  the  present  article,  the  field  energy  a  has  the  dimension 
of  a  reciprocal  area. 

For  absolute  zero,  the  two  entropy  derivatives  can  be  given  explicitly  In 
terms  of  trlgonomet^rlc  functions  as  In  Appendix  C.  Our  analyses 
of  the  case  of  absolute  zero  can  be  summarized  as  follows: 

First,  level  broadening  eases  the  abrupt  variation  of  the  chemical 
potential  as  expressed  by  the  dimensionless  variable  1^ defined  by 
equation  (3.5)  vdien  I/y^  in  the  same  equation  changes.  The  latter  variable 
represents  either  the  electron  density  or  the  magnetic  field.  This 

easement  Is  understandable  because  In  the  presence  of  broadening,  the 
chemical  potential  does  not  jump  from  one  sharp  Landau  level  to  another. 

Actual  broadening  effects  on  the  chemical  potential  must  be  seen  numerically. 
Second,  level  broadening  brings  back  the  Ideal  linear  specific 
heat  which  vanishes  In  strong  magnetic  fields  If  there  Is  no  broadening. 
Interestingly,  this  Important  change  In  the  electronic  specific  heat  takes 


^ace  irrespective  of  the  magnitude  of  broadening.  If  It  exists.  In  a  strong 
magnetic  field  and  In  the  absence  of  level  broadening,  the  electronic 


specific  heat  varies  exponentially  with  the  temperature.  Since  S^ 

In  the  denominator  of  equation  (3.1)  Is  essentially  the  specific  heat,  such 
a  change  from  the  exponential  to  linear  variations  causes  a  considerable 
reduction  In  the  magnetothermal  effect. 

Let  us  now  examine  broadening  effects  more  explicitly. 


5*  Results  and  Discussions 

j 

In  our  grand  ensesible  approach,  the  two  relevant  entropy  derivatives 
with  respect  to  H  and  T  require  the  evaluation  of  the  corresponding  derivatives 
of  the  chemical  potential.,  For  low  temperatures,  these  are  obtained 
by  differentiating  equation  (3.4). 

The  Important  ratio  which  determines  the  magnetothermal  effect 

depends  on  how  the  chemical  potential  u  is  related  to  the  experimental 

variables  such  as  electron  density  n  and  the  magnetic  field  H.  Note  that 

equation  (3.4)  represents  the  chemical  potential  through  y  ^  ■  w/WgH  as 

a  function  of  the  ratio  2iin/p-H  ■  y  Since  their  relation  Is  fundamental 

B  O 

!  ^ 

In  our  theory,  we  have  Illustrated  In  figure  1  the  case  with  T  ■  0.3. 

2 

In  this  graph,  the  three  curves  correspond  to  C  -n  /a  0,  1.5  and  5.0 
respectively.  The  first  curve  represents  the  case  of  absolute  zero  at  which 
the  variation  Is  the  sharpest.  If  there  Is  no  broadening,  the  variation 
at'  0  K  becoaies  zlg-sag  as  in  our  previous  work  (Shlwa  and  Islhars  1983).. 

As  the  temperature  Increases,  the  parameter  C  Increases,  The  case  C  ” 
can  be  considered  Intermediate,  corresponding  roughly  to  1.0  K  and  2  T. 

Such  an  Intermediate  case  is  important  for  experiments.  The  curve  for 
C  *  5.0  Is  almost  straight,  indicating  a  high  temperature  an'd  low  field 
relation. 

.  I 

2 

Figure  2  Illustrates  In  •  the  intermediate  case  of  C  “  w  /o"  1.5,  how 

J 

changes  in  l/y  take  place  about  odd  integral  valuesof  1/y  when  the 
broadening  parameter  ?  Is  changed.  As  we  see,  the  smaller  f"  the  sharper 
the  variation  of  l/y  about  such  points.  The  case  F  ■  3  is  almost  straight 


as  in  the  case  of  C  -  5.0  in  figure  1.  That  is,  both  T  and  C  play  similar 
roles  in  bringing  l/y  close  to  J-/Yq* 

Figure  3  shows  the  magnetothermal  effect  at  two  temperatures  corresponding 

2 

to  IT  /o»  5.0  (left  ordinate)  and  *  1.5  (right  ordinate)  but  at  the  same 

broadening  parameter  ^  ■  0.3.  The  arrows  indicate  the  ordinate  to  be  used. 

These  two  curves  should  be  compared  with  flgi^res  3  and  4  in  I.  We  learn 

2 

that  the  effect  of  broadening  is  strong  in  the  low  temperature  case  of  ir  /o“ 

1.5  where  the  amplitude  is  reduced  nearly  1/100,  but  is  not  very  significant 

2 

in  the  high  temperature  case  of  rr  /a*  5.0.  The  former  corresponds  roughly 

to  1  K  and  2T.  The  region  of  the  abscissa  corresponds  electron  density 
12  -2 

of  order  10  cm  at  2T.  Also,  comparing  the  two  curves  in  figure  3 
with  each  other,  we  learn  that  the  oscillations  are  more  sinusoidal  in  the 

higher  temperature  case. 

/ 

In  view  of  the  coupling  between  the  temperature  and  broadening  effects, 
we  have  Illustrated  in  figure  4  two  oscillations  for  absolute  zero.  The 
left  ordinate  corresponds  to  the  case  ?  ■  0.3  while  the  right  to  3.0. 

We  find  that  the  amplitude  is  reduced  roughly  by  a  factor  of  1/30  due  to 
the  Increase  in  broadening.  Associated  with  the  amplitude  reduction,  the 
oscillating  pattern  is  also  changed  in,  an  interesting  way. 

It  is  important  to  observe  in  the  above  two  figures  that  the  nodes  appear 
close  to,  especially  odd.  Integral  values  of  1/Yq  '  at  low  tempera¬ 

tures  or  small  and  that  the  period  stays  constant: 

A(1/y^)  -  2  (5.1) 

regardless^broadenlng.  On  the  other  hand,  the  pattern  of  the  oscillations 
depends  on  broadening.  The  changes  about  odd  Integral  values  of  the 
abscissa  are  sharper  than  those  at  even  values.  However,  note  that  the 
phase  depends  on  the  g-factor.  We  have  chosen  g  “  2  for  theoretical  con¬ 
venience,  but  the  case  g  *  0  is  also  interesting.  We  shall  comment  on 
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the  phase  In  more  detail  shortly. 

The  oscillations  become  more  sinusoidal  and  the  nodes  deviate  from 
(odd)  Integral  values  when  the  temperature  Is  high  or  the  field  is  low, 
as  in  the  case  without  broadening  which  was  discussed  In  I.  A  single 
parameter  a-  a  /kl  plays  an  Important  role  In  this  respect. 

‘  It  Is  Important  to  remark  that  the  presence  of  broadening  causes 
the  zero  temperature  limit  a:-»-  •>  different  from  the  case  without 
broadening.  We  have  already  pointed  out  that  broadening  restores 
the  linear  specific  heat.  Since  In  the  absence  of  broadening,  the 
specific  heat  Is  exponentially  small,  an  axparlWitttal  determination  of  the 
'  temperatute  variation  of  the  specific  heat  becomes  very  Important. 

\  Since  the  specific  heat  depends  sensitively  on  broadening,  we  have 

Investigated  the  zero  temperature  limit  of  the  magnetothermal  effect. 

Our  results  for  ■>  0.3  and  3.0  are  illustrated  In  figure  4.  The 

I 

arrow  Indicates  the  ordinate  to  Jbe  used.  'in  the  former  cese^r  steep 
changes  occur  near  odd  Integral  values  of  1/Yq.  Its  oscillating  pattern 
Is  significantly  distorted  from  the  sinusoidal  type  which  Is  observed  for.' 

■  3.0.  Hence,  by  observing' such  pattern  changes,  Ijfvel  broadening 
may  be  assessec'.  , 

As  In  figure  3,  the* period  of  oscillation  stays  constant.  Together 
with  the  position  of  the  nodes  at  low  temperature  and  brodaening,  the 
effective  Bohr  magneton  may  be  determined  from  the  magnetothermal 
oscillations.  ' 

However,  for  actual  2D  electron  systems  it  is  necessary  to  extend 
our  consldei^atlon.  First,  let  us  examine  the  choice  of  g  -  2. 

This  Is  a  special  case  which  simplifies  analytical  expressions,  but  in  actual 
Inversion  layers,  the  effective  g  factor  is  enhanced. 

In  general,  the  spin  factor  enters  the  grand  partition  function  as  a 

I 

'  phase  factor,  and  tlic  arbitrary  g  case  can  bo  generated  from  the  case  of  g  ■  0 
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in  accordance  with 


costt/y  -»  4  (co8tii(^  +•§-)+  costii(-  -  -I”)} 

2  Y  ^  Y  2 

On  Che  other  hand,  the  case  g  »  0  is  Biiallar  to  the  Idealized  case  except  Chat 
the  phase  in  terms  of  the  variable  l/y^  is  shifted.  For  example,  in 

the  case  corresponding  to  nig^  X  the  sharp  increases  will  take  place  for  g  ••  0 
at  even  integral  I/y^.  Other  than  that,  there  is  no  significant  change  in 
Che  oscillating  pattern.  This  is  understandable  because  the  oscillations 
‘are  primarily  due  to  the  electron's  orbital  motion. 

The  electrons  in  Si  inversion  layers  have  the  effective  mass  m  -  0.19^^. 
The  affective  Lande's  factor  g  and  the  effective  mass  modify  the  phase 

Z'' 

factor  in  a  combined  form  of  g  'm/2ra^.  This  results  in  a  phase  shift  of 

t  0.3  approximately  if  g  is  3  and  Che  effective  mass  is  0.19ffl  . 

o 


A  comnent  must  be  given  on  the  replacement  of  the  6  functions  by 
Lorentzian  functions  in  view  of  the  work  of  Ando  and  Uemura  \197A)  in  which 
an  elliptic  form  was  adopted.  They  arrived  at  this  form  for  the^  central 

part  of  each  Landau  level  by  neglecting  Landau  level  couplings.  On  the 

. 

other  hand,  Gerhardts  (1976)  arrived  at  a  Gaussin  form  for  level  broadening, 

.  A 

In  what  follows  in  the  present  article,  we  remark  on  the  case  of  elliptic 
broadening. 

-  With  a  proper  normalization,  the  use  of  an  elliptic  form  amounts  in 


our  present  theory  to  hav^*'^ 


o 


% 

vherc  t  1«  generally  different  from  t  of  the  Lorentzlan  case. 

Adopting  the  above  expression  snd  following  the  sane  procedure  as  before 
we  find  that  the  c  dependence  of  not  changed  at  all  but  the 

danplng  factor  of  the  Lorencslan  case: 


Is  replaced  by 


Alt 


?’k 


where  Is  the  Bessel  function  of  the  first  kind,  and 

"If/TS? 

Since  ?  and  are  different,  a  direct  comparison  of  the  two  cases  Is 
sotsewhat  difficult.  However,  as  we  remarked  before,  the  scattering  time 

..  'V 

■  t/2.  Hence,  It  la  meaningful  to  choose  t  «  t/2.  Then,  both 
elliptic  and  Lorentzlan  forms  will  have  the  same  height.  ‘  Figure 
S  Illustrates  the  magnetothernal  oscillations  for  the  elliptic  form 
•  0,6  and  •  1.5. 

While  the  parameter  choice  u  /a*  1,5  for  this  elliptic  case 
Is  the  same  as  the  corresponding  one  In  figure  3,  it  is  Interesting  to 
observe  that  the  oscillating  pattern  Is  similar  to  the  curve  for  P'  ■  0.3 
In  figure  A.  The  nodes  and  period  appear  as  in  the  Lorentzlan  case, 
but  the  amplitude  Is  much  larger  in  tfie  present  elliptic  case  due;  to  the 
lack  of  a  long  tall  In  the  elliptic  broadening  function. 

Although  It  Is  somewhat  beyond  the  scope  of  the  present  article  to 

broadening 

discuss  In  dctsll  the  sbovc  two  cases,  both  cases  seem  to  yield  similar 

A 

results  If  the  respective  broadening  parameter  Is  chosen  suitably. 


In  order  to  apply  the  present  theory  to  actual  systems,  a  question 


reawlns  as  to  the  observability  of  the  effect.  Although  discussions  of 
cxpcrlsMntal  problens  are  also  beyond  the  scope  of  the  present  article, 
we  try  to  conclude  this  article  by  addressing  this  question..  Note  in  this 
resepct,  three-dimensional  cases  have  been  observed. 

For  this  purpose,  let  us  consider  a  model  MOSFET  (metal-oxlde-fleld 
effect  transistor),  consisting  of  a  50A  thick  Ti  gate,  a  4,000  A  thick 
SlOj  layer,  a  50  A  thick  Inversion  layer  and  a  500  um  Si  layer.  The 
surface  areas  of  all  these  components  are  assumed  to  be  the  same,  because  then 
the  area  does  not  affect  our  following  estimate.  Hence,  let  .us  take  a 
unit  area  for  all  the  components. 

Let  us  consider  the  case  in  which  the  electron  temperature  is  dropped  by  a 
few  ten  percent  due  to  a  field  change  of  a  few  percent  at  0.1  K. 

Under  the  adiabatic  condition,  the  energy  thus  released  by  the  electrons 

Is  expected  to  be  transfered  into  the  system.  For  a  typical  electron 

12  -2  about  2 

density  n  ^  10  cm  ,  the  heat  energy  thus  generated  is  3x10  cal/cn  . 

On  the  other  hand,  the  heat  capacities  of  all  the  components  are 
estlSMted  for  0.1  K  based  on  the  and  and  T  laws  as  follows: 

-  3.9x10“^^  cal/K 

Si02-  S.lxlO"^^  cal/K 

C  ,  -  1.5x10”^^  cal/K 

el 

C..  -  6.9x10"^^  cal/K 

dX 

The  inelastic  scattering  time  is  of  order  10  *  sec.  sn  cliat  we  assume  equilibrium 
established  rather  quickly  under  the  adiabatic  condition. 

Prom  the  above  figures,  the  heat  capacity  of  the  MOSFET  is  estimated  to  be 


Ixl0~  cal/K  at  6. IK.  If  the  above  heat  energy  Is  entirely  absorved  by  the 
systea*  the  temperature  change  can  be  around  0.03  K.  Since  we  have  neglected 
the  specific  heat  of  the  thermometer  and  its  leads,  the  actual  temperature 
change  can  be  smaller.  On  the  other  hand,  by  reducing  the  temperature 
and  Increasing  the  magnetic  field  and  its  change,  one  can  expect  larger  temperature 
changes.  Therefore,  the  SMgnetotherasl  oscillations  may  be  observable. 

At  least,  two-dlmensionsl  systems  are  more  favorable  than  the  case  of  bulk 
as  far  as  the  energy  conversion  is  concerned. 
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Appendix  A  Evaluation  of  0^ 


The  Poisson  summation  formula  states: 


?  f(n)  *  I  dsf(8)  +  2  I  I  < 


dsf(s)  cos(27rks) 


(A.l) 


The  series  in  equation  (2.10)  Is  then  expressed  by 


.  ^  (c)  ■  ~~  {■rr+2n  f  ^-2'ffrk 


2n 


k-1 


(A.  2) 


where  a^  and  T  has  been  defined  by  equation  (2.13).  This  dimensionless 

D 


boradenlng  parameter  Is  related  to  the  Dingle  temperature  T^  by 


Fw 


*0  kll 


(a;'3) 


U  being  the  cyclotron  frequency, 
c 


Equation  (2.8)  yields 

iSjd.Jal  r  Jblll  .  2sl  (*•«) 


k-1  k^ 


ir^  k-1  k^ 


Introducing  this  expression  into  equation  (2.7)  and  writing 


O  O  '*' 


(A.5) 


we  obtain 


f”  2  4  3f 

o  ' 


(A.  6) 


where  D(r)  has  been  defined  by  equation  (2.13).  Note  that  the  series  In 

>  ■ 

equation  (A.4)  are  characterized  by  the  reduction  factor  W  defined  by  equation 
(2.12).  The  functions  D(r)  and  W(r)  are  related  with  each  other  by 


D(r) 


-- 

J 


w(r) 


dt 


tn(H-t)  V 


(A.7) 


We  also  find 


vs; 


%f«  9«t 


{■it,  .  I'd.  .  LU.liy^r(i-ni 


Tt6 

sinh(Ti6) 


(A. 13) 


H«nc«, 


.’'It. 

co,(-)  ^2,^ 

^  “  So  2 

.inh(2Ji, 


(A.15) 


Therefor** 


fl  (W) 
o 


.  ^  cos<i!k/Y) 


nB  k-1 


k  sinh(ii  k/a) 


(A. 16) 


The  eppeerenc*  of  the  exponential  factor  exp(-2iirk)  is  charact’erlatic  of 
Lorentalan  broadening.  Note  that  thia  factor  decreases  rapidly  witn  k. 
Adding  the  two  contributions,  we  arrive  at  a  basic  formula: 

e 

n  (u)  .  -  li  -  ^  f 

°  3B  B  £“l  I  sinh(Tr^f/a) 

^cordlngly,  we  can  construct  the  grand  potential  for  arbitrary  g 

following  equation  (2.2). 


(dJ7) 


•  O.  •  ^  s  *  » 


tr. 


Appendix  B  Coiiparliion  wlttl  Dingle’s  Original  Method 


(Dingle's  original  method  amounts  tc  replacing  equation  (2.6)  by 


flP(U) 

o 


n-0 


ln[Ue^-^°^] 

Ce-(n4 


(B.l) 


where  the  suffix  D  is  attached  for  distinction.  Note  that  the  lower  bound 
is  Dingle  gave  a  physical  ground  for  this,  although  not  very  convincing. 

The  n*-sum  can  be  replaced  by  an  integral  based  on  Poisson’s  sum  rule.  The 
result  is 


_  r*  /*  asiks 

(Piv)  I  M*  — - 5—2  ln(l+e 

o  ■  *  k—  i  f 

••• 


By-aoc^' 


(B.2) 


^  's- integral  can  be  performed,  but  results  in  divergences  in  the  c>integral 
such  as  seen  in 


Neglecting  all  these  divergences,  and  after  a  somewhat  lengthy  calculation, 
we  arrive  at 

nP,„,  ..»%(?  ■ 

•  1-1  ,2, 

t  8tnh(-i^) 


,1  -.-2n  [2(n-l)],l  .n, 

+2a[g,(— )+  £  o  g„  (— )]} 


’2'2Y'  n 


(B,3) 


where 


gjl*)  "  I*  j 


-•  •  .•  -•  .•  _■  -■  _•  •  .  .  V.  “  •  •  . 


; .  .1  . 


.  • .  •  . '  .*• 


z  i  -  L.- 

P  3x  IT  ^2^J,2 


(1-2^'^**)  B 

bp  -  <->"  •  (W" 


B  -  <-1l^)"  b 
n  n 


■  ’o'’''^’  -  T<ll>\'‘‘'T.  - 


2Y  «  o - *  “ 

and  B  represents  the  Bernoulli  number  such  as  B  ”  -1,  B,*  1/6,  etc.  Note 

o  o  1 

that  the  first  term  1/2  of  g  (x),  when  Introduced  Into  the  double  Integral 

o 

for  g2(x),  still  causes  a  divergence.  We  must  neglect  this  divergence  also. 
Then,  for  the  case  of  ut>>1,  the  nonoscillating  part  of  equation  (B.3)  yields 
the  sane  diamagnetic  terms  as  Dingle's  (See  his  equation  (4.2)).  Also, 
we  can  show  that  (B.3)  is  reduced  to  the  correct  expression  in  the  limit 


In  our  present  somewhat  Improved  method,  the  replacement  of  Dirac's  5 
function^  by  l^rcntzlnn  funcltons  Is  made  after  '2^(ii)  Is  integrated  twice 
by  parts  as  In  equation  (2.7)  where  the  lower  bound  of  the  c  integration 
Is  kept  to  be  (H-.  We  then  obtain 


4Aa^  r  2ac-Ti  f'l.k  -27irlkJ 

it6  (e  +1) 


r-  f 

o+  0+ 


_  ^  2TrikG  "I 

dc  e  J 


where  means  summation  excluding  the  term  k  =  0,  and 


g.  (e) 


( «'  r 


de"'g  (c") 


Under  the  dHvA  condition  that  y<<1.  8^(1/2y)  is  a  power  series  in 

Y^.  If  1/Y»r,  i.e.,pT>>l  then 

g^(c)'''  1 
o 

R^(c)'^  C^/2 

in  the  region  of  c  inunediately  close  to  1/y  .  Hence,  we  arrive  at 
equations  (A.9)!  and  (A.16).  Note  in  equation  (A.  9)  ,  D(r )  appears  in 
the  nonoscillating  part,  reducing  the  magnitude  of  the  diamagnetic 
susceptibility. 


Appendix  C  Treatment  for  Absolute  Zero 

In  this  section,  we  consider  the  limits  a-»«>  and  B/T-^  and  try  to 
derive  explicit  limiting  formula.  For  this  purpose,  we  introduce 
simplifying  notations  such  as  '  ’ 


?  «  .  .  "Ir  aslnx  i  . 

r  —  sin  nx  -  tan  — 3  -  S_^(x;a) 

n«l 


?  a*^  sin  nx 
n*l 


asinx 

l-2acosx+a 


-  *  S  (x;a) 
2  O 


(C.l) 


»  n 

£  a  cos  nx 
n«l 


acosx-a 


l-2acosx+a 


—  -  C  (x;a) 
2  o 


9  n  .  a[  (1+a  )cosx-2a]  ^  . 

I  an  cos  nX  *  - -  •*  C.  (x;a) 

n*l  (l-2acosx+a^) 


T  a"n  sin  nx 
n»l 


a(l-a^)sinx  ..  , 
- r-r  «=  S,(x;a) 

(l-2acosx+a  ) ^ 


for  a  <  1. 


Hence,  Y  versus  Y  relation  (3.7)  reads  as 
o 


~  -  s  ,  (-;W)  ,  (a-H®) 

Y  Y  ’f  -1  Y 
'o 


(C.2) 


It  is  easy  to  show  that  equation  (3.13)  leads  to 
2 

^  {fCj(~W)+S^(~W)ty  +2S_j(^;W)+2fs^(~W)][l+2C^(~W)]“^},  (a-K»,6/T-»<») 

(C.3) 

We  can  also  show  that 


(C.4) 


We -remark  that  this  entropy  derivative  yields  the  specific  heat  as  follows: 


27 


2  2 

c  -  S  -  t{1+2C  Aw)}  ,  (per  area)  (C.5) 

A#H  Ti  T  O  Of 

Note  tlxat  the  first  ten  is  the  correct  linear  specific  heat  in  the 
absence  of  nagnetic  field.  This  ten  is  canceled  out  if  there  is  no 
broadening  of  the  Landau  levels.  In  fact,  the  specific  heat  decreases 
exponentially. ■  The  reappearance  of  the  linear  specific  heat  is  due  to 

I 

the  introduction  of  level  broadening  which  in  effect  suppresses  the 
oscillating  tens.  ^  The  reappearance  of  the  linear  specific  heat  is  in 
accord  with  Isihara  and  Kojima  who  treated  the  oscillating  terms  to  be 
small  even  without  broadening.  However,  their  assumption  is  valid  only 
for  relatively  high  temperature  and  low  magnetic  field. 

These  limiting  formulae  can  further  be  simplified  if  l/y  is  an 
integer:  It  U  jt  1  and  l/y  -  m  ,  an  integer,  we  find 

S  ■  0  n  ■  -1,  0,  1 

n 


c  .  _ _ !-i>-  , 

*  —  ^ 

Hence, 

”  ^  [i-(-)"w]^ 


S 


T 


3a 


1- (-)"w 


We  obtain 


2wr  ^ 

l-K*  " 


/even 
“  'odd 


(C.6) 


(C.7^ 


(C.R), 


(C.9) 


Although  those  expressions  Involve  series,  their  convergence  is  very 
good  in  general  due  to  the  presence  of  the  convergence  factor 
appearing  through  and 


i 

i 


References 


Ando  T  and  Uemura  Y  1974  J.  Phys.  Soc.  Jpn  ^  959 
Brailsford  A  D  1966  Phys.  Rev.  139  456 

/m 

Dingle  R  B  Proc.  Roy.  Soc.  (London)  A  211  517 

/>■ 

Gerhardts  R  R  1976  Surf.  Scl.  ^  227 

Islhara  A  and  Kojlma  D  Y  1979  Phys.  Rev.  B  1£  846 

_  and  Shiwa  Y  1983  Solid  State  Coramun.  48  1081 


_  1984  J.  Phys.  C:  Solid  State  Phys;  17  5075 

_ _  1984a  Proc.  17th  Interflat.  Conf.  on  Low  Temp.  Phys 

ed.  by  Eckem  U  et  al  (North-Holland.  Amsterdam)  j)«,..fi02^  j’ 7  / 

_ _  ,  Tsai,  J  T  and  Wadati  M  1971  Phys.  Rev.  A  3  990 

McConbe  B  D  and  Seidel  G  1967  Phys.  Rev.  ^  633 
Shiwa  Y  and  Isihara  A  1983  Phys.  Rev.  B  4743 
_  1984  Solid  State  Connun.  49  497 


Zawadski  W  and  Lassnig  1984  Surf.  Sci.  142  225 


Figure  Captions 


Fig.  1. 


Fig.  2. 


Fig.  3. 

/ 


Fig.  4. 


y  Fig.  5. 


Variation  of  l/y*  as  a  function  of  l/v  ■  2Ttn/iJ  H 

(H  “  1,  2in  ■  1)  at  ?  "  2irr“  0.3  for  three  different  temperatures. 

Here,  u  is  the  chemical  potential,  n  is  the  electron  density  and 

lij  la  the  effective  Bohr  magneton,  r  being  a  broadening  parameter. 

Effect  of  broadening  on  1/y  which  Is  the  chemical  potential  In  the 

2 

units  of  the  field  energy  at  constant  n  /a-  1.5,  where  a  -p  H/kT. 

B 

This  case  corresponds  approximately  to  1  K  and  2T. 

Magnetothermal  oscillations  for  a  fixed  broadening  parameter 

^  ■■2nr-  0.3  at  different  temperatures  corresponding  to  5.0 

2  ^ 

(left  ordinate)  and  ir  /a-  073  (right  ordinate).  The  latter 

represents  roughly  the  case  of  2T  and  IK  In  the  region  of  electron 
12  2 

density  of  order  10  /cm  . 

Broadening  effects  at  absolute  zero.  Left  ordinate:  ?*■  0.3. 
Right  ordinate:  3.0.  The  arrow  In  each  curve  represents  the 

ordinate  to  be  used. 

m 

Magnetotheral  oscillations  for  an  elliptic  density  of  states 

with  ?'  -  2nr*“  0.6  and  for  w^/a  -  1.5.  T'  Is  the  broadening 
parameter  In  this  case. 
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